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ABSTRACT 


An aeroelastic investigation of horizontal axis wind turbines 
is described. The study is divided into two simpler areas, namely, 
the aeroelastic stability of a single blade on a rigid tower, and 
the mechanical vibrations of the rotor system on a flexible tower. 
Some resulting instabilities and forced vibration behavior are 
described. 


INTRODUCTION 


The aeroelastic analysis of horizontal axis wind turbines may 
be divided into two convenient areas, namely, (a) the investigation 
of the aeroelastic stability and response of a single blade on a 
rigid tower, and (b) the investigation of the mechanical stability 
and vibrations of the rotor system on a flexible tower. With a 
reasonable understanding of the behavior in these two areas, the 
more complex problem involving a completely coupled blade-tower 
aeroelastic system can be better understood. 

The present report deals with some work performed for ERDA in 
the two areas cited above. The aim has been to investigate the 
basic mechanisms and to reduce these analyses to forms in which 
they can be readily used for trade-off analyses. 

AEROELASTIC STABILITY OF SINGLE BLADE 


Using the equivalent hinge concept, the nonlinear equations of 
motion were written for large deflections of a rigid blade, with 
elastic hinges at the root, rotating about a rigid tower. The blade 
was assumed to have a flapping (3) / lagging ((j>) , and feathering (6) 
degree of freedom, and the effects of preconing of the blade rela- 
tive to the feathering hinge (5), and of the feathering hinge rela- 
tive to the plane of rotation (Sjj) / were included. Quasi-steady 
aerodynamic forces were assumed. The equations were first solved 
to obtain the static, steady state displacement in flap (3q) and 
lag ((Jjq) . Then the equations of motion were linearized by assuming 
small perturbations about the steady state flap and lag displace-- 
ments found above. The resulting linear equations take the form 

[A]{x} + [dHx> + [kHx} = 0 (1) 
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where {x} represents the perturbation displacements in 3 » <j)/ 6 
[A], [D] , [K] represent the 3x3 mass, damping, and stiffness matrices 

respectively. These equations were investigated for stability by 
recasting them as standard first order equations, and extracting 
the eigenvalues p^ = + i bn using standard procedures. The ^ 

resulting a^'s indicate the damping characteristics, and the b^^ s 
the frequencies of the resulting oscillations. 

A wide range of parametric variations are being investigated. 
Figure 1 shows some typical results from this analysis showing 
effects of preconing on the stability margins. The results 
presented in the form of modal critical damping ratio ^ ~ — ^n' n 
since it is believed this more clearly indicates the strength of the 
instability (or stability) than simple stability boundaries. It is 
seen here that positive droop angles (6) can particularly decrease 
the stability margin. 

The analyses described above have been based on equivalent 
hinges in order to provide a rapid and reasonably accurate tool for 
predicting stability characteristics. However, a more elaborate 
program is being developed which allows for distributed flexibility 
of all three hinges, and this program will be used to validate the 
results using equivalent hinges. 

The linear aeroelastic stability analyses have been supplement- 
ed by some nonlinear, large deflection analyses again using equiva- 
lent hinges, to determine the nonlinear limit cycle behavior of 
these instabilities. The harmonic balance method was applied and 
the resulting nonlinear equations were solved by the Newton- Raphson 
technique. Both the self-excited aeroelastic instabilities and the 
dynamic response to gravity and wind shear loads were investigated. 
Figure 2 shows the nonlinear limit cycle behavior of a predominantly 
flapping- torsion type instability as rotation speed is increased. 
The nonlinearity tends to be of a softening type, i.e., steady limit 
cycles can occur somewhat below the linear critical rotation speed 
if a large enough disturbance is given to the system. 

MECHANICAL VIBRATIONS OF ROTOR- TOWER SYSTEM 


The mechanical instabilities and vibrations that may result 
from the interaction of the flexible rotating blades with the base 
motions of the supporting tower were also investigated. The sim- 
plest of these interactions involves the coupling of the tower side 
motions (qL) with the blade lagging motions (<}>), which can give 
rise to a strong mechanical instability (the "gound resonance effect 
of helicopters) , as well as significant forced vibration effects 
due to blade unbalance and gravity forces. Using the equivalent 
hinge concept, the equations of motion for small vibrations of such 
a rotor- tower system are 

i ^2 2 i 

M,q^ + c^q^ + k^q^ + I ^ 1 S.sin 

Xj Li L Li L L X 1 X 

cos4j^ + 1^$^ + = g[S^sini|;^ + S^c()^cosijj^] 


i = 1,2, .. .N 


( 2 ) 
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where ijj. = Jit + 27r(i-l)/N represents the angular position of the 
i"*^^ blaae, g is the acceleration of gravity, is the equivalent 
tower mass of the tower and rotor, and S- and, represent the 
static moment and moment of inertia of the i^” olade about the lag 
hinge. 

The mechanical stability of this rotor-tower system can be 
found by neglecting gravity and blade unbalance forces 
(g = 0 , Si = const), and introducing the multiblade coordinate 
transformation 


(j)^ = a(t) sim|;^ + b(t) cosij;^ 


(3) 


This eliminates all the periodic coefficients in Eqs. (2) for N>3, 
and the stability of the resulting equations in q^, a, b, can be 
investigated in the standard manner. Figure 3 shows the results 
for a typical 3-bladed rotor with varying amounts of critical 
damping ratio in the tower and in the blades. In general adding 
damping stabilizes the system, but under certain circumstances, the 
addition of unequal amounts of damping can destabilize the system 
from the no damping case. For a 2-bladed rotor, effects are similar, 
but with additional small instability regions. 

Effects of small, static unbalance can be found by assuming 
that Sj = S + AS, S 2 = S 3 = ... = S. The right hand side of 
Eqs. (2) then reduces to the forcing function 9?- S sin J^t. By neg- 
lecting the small AS variations in the periodic coefficients, and 
introducing again the multiblade coordinate transformation, Eq. (3) , 
the basic equations reduce to the forced harmonic response of a 
conventional spring mass system involving the coordinates qj^, a, b. 
The resulting tower and rotor responses become 




qj^(t) = qj^g sinJ 2 t + qj^^ cosQt 
b +a b +a 

- sin 2 J 2 t + ( 2 ) cos 2 fit 


(4) 


where ag, a^, bg, b^, are the sine and cosine responses of the coor- 
dinates, a, b. Figure 4 shows the amplitude of the blade lagging 
response for a typical case. It is seen that strong 2/rev 
resonance responses occur in the blade when the rotation speed 
equals 1/2 the lagging natural frequency or when Q, equals the 
tower natural frequency . Adding damping to the blade or to the 
tower reduces the respective amplitude accordingly. 

The effects of gravity loading can similarly be investigated 
by examining Eqs. (2). Generally, gravity acts directly on each 
blade without much tower interaction. Strong 1/rev resonances may 
occur when J2 equals the lagging frequency w. and weak .5/rev para- 
metric resonances may possibly occur when Ji^equals twice 

More comprehensive rotor-tower interactions were investigated 
by including the tower yawing motions (0^) and the blade flapping 
motions (3) in the previous analyses. The mechanical instability 
regions occurred at slightly lower rotation speeds 
than before, but the general mechanisms were not changed signifi- 
cantly. These additional degrees of freedom 
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also affected somewhat the static unbalance loads. Of more signi- 
ficance was the direct effect of wind shear and tower shadow on the 
blade flapping loads when these additional motions were introduced. 
The investigation of these mechanical vibrations are continuing. 
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DISCUSSION 

Q. Comment; It is the flywheel resonance and mechanical insta- 
bility problems that have primarily dictated natural fre- 
quencies of the tower well above the operating range of the 
rotor. 

A. Yes. These factors play important roles, together with the 
forced vibrations due to rotor static unbalance, gravity 
loads, and air loads. Generally, I believe, to avoid signi- 
ficant "ground resonance" type mechanical instability, the 
blade lagging frequency (in-plane vibrations) would need to 
be above the rotation speed of the rotor. Low natural 
frequencies of the tower usually give rise to large forced 
vibration as the rotor speed passes through the tower natural 
frequency (and also through one half the tower natural 
frequency for two bladed rotors) . 
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Figure 1 - Effect of precone and droop on aeroelastic stability of lag mode. 


57 




Flutter solution (nonlinear limit cycles). 
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Rotation speej> ratio ~ 

Figure 4. - Forced response of rotor to static unbalance. 
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